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■ Abstract 

(N 

We model the classical transmission of a massless scalar field from a source to 
i a detector on a background causal set. The predictions do not differ significantly 

■ from those of the continuum. Thus, introducing an intrinsic inexactitude to lengths 
and durations - or more specifically, replacing the Lorentzian manifold with an 
underlying discrete structure - need not disrupt the usual dynamics of propagation. 

cr 

£jy 1 Introduction 



Despite the variety of approaches the the problem of quantum gravity, it is hoped by some 
that agreement can be forged on suitably generic consequences of the as yet unknown 
theory. This would offer hope of deriving generic predictions of quantum gravity in a 
| more or less heuristic fashion. However, deriving such predictions often turns out to 

demand a greater level of specificity, whereby disagreements return. For instance, a 
majority of workers would probably agree that the diffcrcntiablc manifold structure of 
, spacetime will break down near the Planck scale to be replaced by something of a more 

discrete, "quantised" or foamy nature. When potentially observable consequences of these 
general ideas are sought, however, the consensus evaporates. For example, disagreement 
arises over whether the expected break-down of General Relativity at Planck scales would 
give rise to modified dispersion relations or other Lorentz symmetry violating phenomena. 
Spacetime discreteness is often cited as motivation to consider Lorentz symmetry violation 
(see e.g. pQ), perhaps due to consideration of altered dispersion in simple lattice models. 
There are also quantum gravity inspired models that draw the same conclusion [3J [3] . 
Meanwhile, there are models of discrete spacetime which respect Lorentz invariance [H [51 
[3] and so do not result in Lorentz symmetry violating phenomenology. Another example 
of controversy is the arguments in references [7, 8, 9 which dispute the generic claim of 
reference [101 that any quantum gravitational "fuzzing" of the metric would disturb the 
coherence of light from distant sources by an amount greater than that allowed by the 
observations. 

The lesson is that it is important to test expectations of phenomenological effects 
of quantum gravity against different concrete models in order to determine if they are 
generic. One might indeed be tempted to claim that any "fuzzing" of spacetime properties 
- be it quantum uncertainty in time and distance measurements, lack of definition of short 
distances due to discreteness, or from some other source - would naturally lead to some loss 
of coherence of light from distant sources (for a model inspired by similar considerations, 
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see [TT]). One purpose of this paper is to test this expectation in a definite model. More 
specifically, the model examines this question in the discrete context of causal set theory. 
We present a model in which spacetime metric relations do indeed have an approximate 
character that breaks down at around the Planck scale, but in which no significant loss 
of phase information results. 

One crucial point is that, whereas other models introduce Lorentz violating assump- 
tions, the causal set model allows for Lorentz invariance in the approximating continuum. 
That the transmission of light in this context must be consistent with Lorentz symmetry 
follows already from the results of [5j [6] , and in agreement with this conclusion, the model 
we study below will exhibit the usual, strict proportionality of frequency to wavelength. 
However, we know of no equally general reasons that would ensure the coherence of light 
from distant sources; and so one might wonder whether the underlying discreteness would 
necessarily disrupt the coherence of propagated waves, or wrinkle the wave fronts or alter 
the wavelength, and hence the frequency of the signal. 

A second purpose of this paper is to learn something about other possible effects of an 
underlying atomicity that one might expect to encounter in analogy with more familiar 
examples like the propagation of light through air or other material media. Assuming this 
analogy is valid, effects like scattering and extinction will be present at some level, and 
the question becomes whether one can expect them to rise to the level of observability 
with current technology. 

In this paper we formulate a rudimentary modef] with whose aid one can study some 
of these putative effects, assuming that the microscopic structure which replaces the 
continuum at the Planck scale is the causal set. The specific process we will analyze is 
the transmission of a signal from a source to a detector by a (classical) massless scalar 
field in Minkowski spacetime M 4 . Ideally, one would analyze this process with respect 
to a complete theory of electromagnetism formulated within the context of a complete 
theory of quantum gravity. In practise, however, one can hope that a simple model 
incorporating the elements of discreteness and wave-propagation can illuminate the range 
of possibilities to be expected. Here we study the simplest such model we could devise, 
based on a classical scalar field <fi of zero mass. We do not formulate an independent 
dynamics for <j). Rather we describe its transmission from source to detector in terms of 
the discrete analog for the causal set of the retarded Green's function (a delta function 
on the future light cone) that yields the Lienard-Wiechert potentials in the continuum. 
In the future, we hope to compare the predictions of this model with those of a more 
complete theory of wave propagation in a causal set (for progress in this direction, see 

[HIHQ12CL6]). 

2 A continuum model 

As illustrated in figure [IJ we will imagine an idealised setup consisting — in a continuum 
description — of an oscillating point source of scalar charge q, together with a detector 
of rectangular shape, facing the source and at rest with respect to it in an ambient flat 
spacetimqj . In a system of Cartesian coordinates in which the detector is at rest, let the 
source oscillate about the origin, tracing out the trajectory 

x 1 = acosusx , x 2 = x 3 = , (1) 

1 Some features of the model were previously reported in the proceedings of a conference |12| . 

2 For a slightly more realistic model (at least for the comparison to propagation of electromagnetic 
radiation) the source should be neutral. To that end we will later be supplementing this point-charge 
with a static companion of opposite sign; this cancels the constant, "DC" component of the far field but 
leaves the AC component unaffected. We could also include more worldlines to build up an extended 
source, say in the dipole approximation. On the other hand, we could (since scalar charges need not be 
conserved) simplify our setup even further by keeping the charge at rest but letting its strength q vary 
with time. This would take us farther from the astrophysically relevant, electromagnetic case, however. 
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where a is a negligible fraction of the distance R to the detector. (On average, the source 
is thus at rest with respect to the detector, another simplification which one could easily 
relax.) Assuming that the detector has a resolution time of T, a two-dimensional area A, 
and a thickness d, this gives us in spacetime a detection region T> of volume TAd, which 
we can take to be the rectangular domain [R, R + T] x [R,R + d] x [0, VA] x [0, VA}. (In 
this paper we adopt Planckian units in which c — h — 8irG = 1.) Let us assume further 
that R greatly exceeds A ~ the wavelength of the emitted radiation, which in turn 
is much greater than any detector dimension; and also that the detector's response rate 
is rapid compared to the period of the source oscillations. We have then 

d, \f~A~ « A « ]?, T < A, and a < R . 

Up to small corrections, the spatial distance between emitter and receiver is thus R. 
We remark that, in order to observe coherence effects, the detector actually must be 
operational for a much longer time than T, namely a time comparable to a; -1 . Finally, 
let us take the "output" F of our detector to be simply the integral of the field, <fi, over 
the detection region V. 0. 



R 



source 



It 




a) space 



b) spacetime 



Figure 1: The Source and Detector. In (a) the spatial layout is shown, while (b) is a spacetime 
diagram (with only one spatial dimension represented) showing the source's worldline P, the 
light-cone of a point x on that worldline, and the detector region in grey. 

Before proceeding to our causal set model of the same setup, let us first review the 
treatment of this situation in the continuum. The retarded Green's function for a massless 
scalar field in (3 + l)-dimensions is a delta- function on the forward light-cone: 



G(y,x) 



Airr 



w^6(\y — x\ 2 ) if y is in the causal future of x 
3 otherwise 

' 6(y°-x°-r), 



(2) 



(3) 



where r is the spatial distance from x to y. In terms of G, the field produced by our 
source is given by 



G(y,x(s))qds , 



(4) 



3 This model could be regarded as fairly realistic for radio wave reception. For optical frequencies, the 
time-resolution of realistic detectors is not better than the period of the wave, and the received signal is 
usually compared to another oscillator in order to diminish the effective value of T. In such cases our 
parameter T must be understood as an effective resolution-time characterizing the detector. 
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where P is the worldline of the source, q is its charge, and s is proper time along pQ 
According to the ansatz we have made, the output of the detector is proportional to the 
integral J v d 4 y <j>(y) , and thus can be expressed directly in terms of the source and the 
Green's function: 



F = q ds d 4 yG(y,x(s)) (5) 
Jp Jv 

The two integrals can be done in either order. Integrating first with respect to s yields 
for (j>(y) 

g ds/dt 

AirR 1 - dx/dt ' 1 ' 

where ds/dt and dx/dt are to be evaluated at the retarded time t — R corresponding to 
the detection time t of interest. (Strictly speaking, this expression is correct only at a 
single point within the detector close to (R, 0,0), but we can ignore this caveat, since 
we have assumed the detector to be small compared to R and T small compared to the 
period of oscillation, In ([6]), the oscillating signal resides in the time-dependent factor 
multiplying l/R, which can be understood as just the Doppler effect in disguise. The 
y-integration is now trivial and furnishes for the detector output 



F « J±±Z A TAd , (7) 
V 1 - v AttR ' v ' 

where v = dx/dt is the component of the source velocity toward the detector at the 
relevant retarded time. In a slow velocity approximation, yl + v /\ — v reduces to 1 + is, 
and the varying, "AC" part of the signal is a multiple of sin(wt)/47ri?. 

In order explicitly to subtract off the DC part of the signal and create a more realistic 
model of electromagnetic radiation, we can add a static negative charge, —q, at the origin 
of spatial coordinates. Then, renaming F above as F + and calling the detector response 
to the negative charge, F_ , we have the total detector response 



r ^--= r - +F - = ^m TAd -m TAd (8) 



which, for low velocities is 



F to t a i = is^TAd (9) 

= -^—TA dauj sm(ujt) . (10) 
47ri? v ' x ' 

The detector response is coherent in space as well as in time, so our model captures the 
coherence that spacetime uncertainty might be expected to disrupt. 

Note that the l/R dependence of the signal agrees with the usual falloff of far field 
electromagnetic radiation from a dipole source. Here the signal is proportional to the ve- 
locity of the source however, whereas the signal in the electromagnetic case is proportional 
to the acceleration. We note that we could make the model more realistic by pretending 
that the "field" <p(y) is the electrostatic potential Aq of an electromagnetic field. The 
relevant spatial gradient of Aq = (f>(y) would arise from differential time delay (detectors 
at different spatial positions see the source at different times along its world line) and 
could give rise to an electric field proportional to the acceleration of the source. 



4 We have normalized G so that n G = — <5, where our metric signature is ( — hH — h). Notice that only 
an infinitesimal portion of P contributes to <j>(y). 

5 The fact that different points in the detector correspond to different values of s could not have been 
ignored if we had integrated first over y, however. 
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3 The discrete model 



Now let us consider a causal set model for the same situation. Before turning to the 
model per se let us recall some kinematical results from causal set theory. (For a more 
detailed introduction see [H[T71[T3].) A causal set C, or causet for short, is a locally finite, 
partially ordered set. When a causal set has an approximation by a continuum spacetime, 
the order-relation of C corresponds to the causal order of spacetime, while the number 
of elements in a subset of C equals (up to fluctuations) the volume of the corresponding 
spacetime region in fundamental units. It is conjectured that this information is all that 
is needed to recover the metric, differential structure, dimension and topology of the 
approximating Lorentzian manifold (at scales large compared to the fundamental scale). 

With a nonzero discreteness scale, the correspondence of any Lorentzian manifold to 
a causet can only be approximate, and it seems to be necessary to define it stochastically. 
Although the causet is taken to be fundamental and the manifold is merely an effective 
description of it on large scales, some means must be found of determining whether a given 
spacetime approximates to a certain causal set. To this end we consider a method called 
sprinkling, which produces a random causet from a given Lorentzian manifold by means of 
a Poisson process [18]. One samples, at random, points from the manifold with a density 
of 1 in some fundamental units. (One might term these fundamental units "causet units" . 
They are the analog of lattice units in lattice gauge theory, and one would naturally expect 
them to agree approximately with Planck units as defined above; for the purposes of this 
paper we will take the agreement to be exact. To determine the true factor relating the 
two sets of units would be to determine the fundamental discreteness scale and would be 
a key step in the development of the theory.) A causet can then be constructed using 
these "sprinkled" points as elements, the order-relation being that induced by the causal 
order of the Lorentzian manifold. A manifold is defined to approximate a causal set when 
the causal set is a typical result of the sprinkling process on that manifold. 

Often, instead of finding an approximating spacetime for a given causal set, we are 
given the approximating spacetime, and wish to find a causal set to which it approximates. 
In that case, we can reverse the argument and construct possible underlying causal sets 
using the sprinkling process. If, in the sprinkling process, a property of the resulting 
causal set holds with probability 1, we say that the property holds for a typical sprinkling. 
Similarly, properties that hold with probabilities sufficiently near to 1 can be considered 
typical. Some care is needed to properly define this notion in all cases. However, without 
fully entering into this discussion, if the only quantity of interest is a function x(C) of 
the causal set, the value for a typical sprinkling will be only a small number of standard 
deviations away from the mean over sprinklings (x) . One can imagine cases in which x is 
an "ensemble average" of many variables pertaining to different regions. 

In a Poisson process, the probability to sprinkle n elements into a spacetime region of 
volume V is given (for any measurable subset) by the Poisson distribution, 

V n e~ v 

P{n) = j— ■ (11) 

Because a Poisson process depends only on the spacetime volume-element, it is Lorentz 
invariant. This invariance is exact for d-dimensional Minkowski spacetime, where the 
theorems on the existence and uniqueness of the Poisson process readily establish its 
invariance under all volume preserving linear maps, and in particular under arbitrary 
Poincare transformations (see e.g. reference [T8]). For a proof of invariance that holds 
even for individual realizations, see [6]. Thus, if the causal set is well-approximated by 
Minkowski space at all, then the discreteness does nothing to pick out a preferred frame 
in it. In a curved spacetime, Lorentz invariance is to be understood in the same, local 
and approximate sense that holds in General Relativity. 

For future reference we define the terms "interval" , "chain" , "path" and "link" . Given 
elements p and q in C, the order-interval between them is I(p, q) = {r G C \ p -< r -< q}, 
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where -< is the fundamental precedence relation defining C. A chain is a linearly ordered 
subset of C and a path is a saturated chain, i.e. one which is maximal within the order- 
interval between any two of its elements. A link is a causal relation p -< q between two 
elements, p and q, that is not implied via transitivity by other relations, i.e. there exists 
no third element r such that p -< r < q. 

We now turn to the causet description of the source-detector system. We will describe 
the discrete counterparts of the spacetime, the source and detector, and the model of 
propagation used above, and determine the detector output in this model. In place of 
4-dimensional Minkowski space M 4 , we will have a typical sprinkling of M , and in place 
of the detector region, the subset, T>, comprising all the elements that were sprinkled into 
the region of spacetime T>. The field <fi will be a function from the elements of the causet 



We will ignore, for now, the static negative charge source and study only the discrete 
analogue of the oscillating positive charge. In analogy with the continuum model of the 
source as a point-charge, we may in the causet identify the source with a path P that 
approximately follows the spacetime worldline P of eqn. (JXJ) . Such a path will be "locally 
geodesic" in the sense of approximating a longest chain between any two sufficiently nearby 
elements of P. We will say that an element of P is associated to a point x of P if it lies 
at a larger spacelike distance from x than any other point in P. The sprinkled points 
making up such a P will be distributed along the corresponding continuum curve with 
an approximate spacing in proper time along the curve that is near to one Planck unit 
[T9l [20] . (The ratio of chain- length in the causet to proper length in the continuum does 
not seem to be known exactly in 4-dimensions. Simulations suggest a value somewhat 
greater than 1.2. This gives the average number of elements per length of P and this, 
along with the independence of the sprinkling process in different regions, turns out to 
be all that is needed for our calculations.) Such a locally geodesic path will exist, and 
the spatial distance in the sprinkling from the curve P to the points corresponding to 
elements of P will typically be a small number in Planck units. 

Both of these definitions make use of regions in the spacetime M that we are to 
sprinkle, rather than being defined "intrinsically" by relations and/or fields on the causal 
set itself. Defining them by such means would be more satisfying physically but would 
complicate matters considerably, and if it led only to a change in the position of the two 
regions on or around the order of the Planck scale, that would not affect the results, as 
will become obvious during the calculation. 

In the continuum, the relevant dynamics of the field were completely described by a 
retarded Green's function G(x, y). To define an analogous "propagator" in the causal set, 
we must discover a discrete replacement for G(x,y) that approximates it well on large 
scales, but relies only on the structure of the causal set, without appealing to any extra 
information from the continuum. This idea has already been explored in [21j , [22) (which 
are concerned primarily with the 2D situation) and [23] . 

The task of finding a discrete propagator is made easier by the simple form of equation 
IJ5J, as we now describe, using the link concept defined above. In a sprinkling of M 4 , the 
future links from any given sprinkled point are unlikely to stretch over long proper times. 
Nevertheless, thanks to the Lorentz invariance of the Poisson process, links from e» will 
connect it with an infinite number of other sprinkled points ( "future nearest neighbors" ) 
spaced out along, and just inside of, the future light-cone. This assertion can be proved 
as follows. 

The mean number of links joining an element at point x € M 4 with elements in a 
region 1Z contained within its future is 



where \I(x,y)\ is the spacetime volume of I(x,y) — J + (x) n J l(y), the causal order- 
interval (or "Alexandrov set") between x and y. This expression can be understood as the 



to R. 




(12) 
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sum over infinitesimal volumes d 4 y in 1Z of the probability that an element is sprinkled 
into d y times the probability that this element is linked to the element at x, i.e. the 
probability that no point was sprinkled into the interval I(x,y). For a Poisson process, 
the former probability is d 4 y itself and the latter is e~^^ x ' v ^, from which (|12[) follows. 
Now in four dimensions, 

\I(x,y)\ = ^r*(x,y) , (13) 

where t(x, y) is the proper time between x and y. From this one sees that an element 
that was sprinkled a large proper time from x is highly unlikely to be linked to it; but 
by the same token, elements hugging the future light-cone of x which are close to it in 
proper time are quite likely to be linked to it. 

Indeed, consider an element eo sprinkled at the point x e M 4 . Choose coordinates 
£, 77, 9, if {rj> 0), in which x is at the origin, £ = 0, and the metric is 

ds 2 = -df + fdrj 2 + f sinh 2 rjdfll (14) 

so that £ is geodesic proper time from the origin. The expected number of future links 
from eo terminating in the region between the light-cone of x and the hyperboloid of 
points at a fixed proper time r from x is given by (|12[) with 1Z defined by < £ < t, that 
is by 

" i>T "I />OC 

47T / <i££ 3 e~^ 4 / dr) sinhr) . (15) 
Jo J Jo 

Since the integral over r\ diverges, the expected number of links is infinite, no matter how 
small t is. 

We introduce the causet function 

whenever e < e' and {e, e'} is a link, 

(16) 

otherwise, 

where e, e' € C are causet elements and k is a normalising constant of order 1, to be decided 
later. In the limit of infinitely dense sprinkling, this function becomes a J-function on the 
forward light-cone of e. This means that L(e', e) is a Lorentz invariant discretisation of 
the continuum Green's function, and can be used to define the propagation of the scalar 
field on the causet. 




3.1 Mean Detector Output 

We will calculate the average and variance of the detector output over sprinklings of M 4 . 
Replacing G(x,y) with L(e',e), the calculation of the detector output is very similar to 
that in the continuum, although now it is a random variable. The double integral of eqn. 
([5]) for the output of the detector F becomes a double sum over all elements in the source 
and detector: 

F = <?£I>(e',e), (17) 

eEP e'ET> 

In other words F is proportional to the number of links from elements in the source 
region to elements in the detector region. The causal set is randomly generated by the 
sprinkling process in M 4 , with the source and detector regions as described above. Using 
the definition of L(e', e) in (ITBl . F can be rewritten as 

F = qn( f a(ds) X (d 4 y)a^),y), (18) 

J P JT>x(s)^y 

where s and x(s) are as in the continuum calculation, and a(ds), x{d A x) <^(x, y) are random 
variables in the sprinkling process: the variable a(ds) — 1 if any point of P with proper 
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time in ds has a sprinkled element of P associated to it and is otherwise, x(d 4 y) = 1 if 
the volume element d y contains a sprinkled point and otherwise, and £(x, y) = 1 if the 
causal interval between x and y (not including the points x and y) is empty of sprinkled 
points and otherwise. From this, we can hnd the mean value: 

{F)=qn( [ (a(d S ) X (d 4 y)ax( S ),y)), (19) 

JP JVx(s)^y 

= q K [ [ (a(ds))( X (d*y))(ax(sU)), (20) 

JP JVx(s)^,y 

= qn [ [ daefiye-M'W'rt. (21) 

JP JVx(s)^y 

As explained above, the average number of points in P associated to P per unit length, 
{aids)) / ds, is some number of order unity, which we absorb into n. The calculation is 
unaffected by the Planck scale deviation of the spatial positions of elements in P from 
P. The second line above says that the three random variables are always uncorrelated: 
firstly x(s) and y range over disjoint regions, so cr(ds) and x(d 4 y) are uncorrelated, and 
secondly points sprinkled at x(s) and y have no effect on £(x(s),y). Let us first calculate 

(4>{y))=qn( dse-I^W'*)!, (22) 

J P x(s)^y 

for y e D. 

(Despite the name this quantity is different to the discrete field, which is a function 
of the causal set elements). It is convenient to define r = y 1 — a; 1 (s) and t — y° — x°(s), 
and use the co-ordinates u = t — r, v = t + r to calculate the integral: 

(4>{y))^qK dse-^ uv \ (23) 

where we have used (|13[) and used <C R to approximate the proper time between x(s) 
and y. Let us define R y — r\ u —o. Using this, 



where 



We also have 



au + 2R y , (24) 

(25) 

1 — V 



ds = \/i - v 2 dx° = -y/adu. (26) 
Inserting these into ([23jl yields 



/•OO 

{cj>{y)) w / du exp[-^(cra 2 + 2R y u) 2 } (27) 

Jo ^ 4 

(28) 

since the integral is dominated by the region in which u -C R, where the integral is 
approximately Gaussian. Note that R y is approximately R, because the amplitude of the 
source a <C R and d <C i?. Inserting this result into eqn. (|2ip we obtain our result for the 
mean output, 

^"\fr IE ^ TAd " {1+ ^ TAd - < 29 > 
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in the low velocity approximation, where we have set k = 1/(2v6tt)- We see that the 
mean output in the discrete model matches the continuum result given above. We 
have approximated the spatial distance from source points to detector points as R in all 
cases when finding <f>(y), but this is basically the same approximation we made in the 
continuum case. 

It is of interest to find the size of the correction to the continuum result. The main 
extra approximation used in the discrete case was the approximation of the integral in 
equation (|27|) as Gaussian. Instead of throwing away all but the Gaussian factor, we can 
include the next most significant factor: 

~ f°° TT 

(<j>{y)) « qnV^ / du exp[--(R y au 3 + R 2 y u 2 )] . (30) 
Jo o 

Solving this integral by symbolic computation, and expanding in powers of 1/R, shows 
a first term equal to the continuum result, and a second term proportionally smaller by 
a factor of order 1/R 2 . This R is a large distance expressed in Planck units, and so the 
extra approximation introduced in the discrete case is insignificant. 

Just as in the continuum case, we cancel the DC component of the response by adding 
a static, negative charge, modelled on the causal set as a path close to the origin. The 
mean of the total response is the sum of the means of the responses to the two sources 
and is 

(Ftot^av-l-TAd. (31) 
4ttR 

3.2 Fluctuations in the detector output 

We have shown above how, in our model of propagation, the average value (over sprin- 
klings) of the detector signal does not significantly differ from the signal in the standard 
case. However, this result must be strengthened. We want to make sure that the result 
does not differ significantly not just for the average, but for a typical sprinkling (moreover 
we want to make sure that the model predicts the correct result for many repeats of the 
experiment). To this end we calculate the variance of the signal. We start by calculating 
the variance when the source consists of a single oscillating positive charge. 

As stated above, the detector output F is proportional to the number of links from 
the source to the detector. If the existence of any link was uncorrelated to the existence 
of any other link, we would expect a variance of qn{F) which is small. Correlations do 
exist, however, intuitively they will not be large. The result of this intuition - that the 
fluctuations, as a proportion of the signal, are small for reasonable values of the parameters 
- is correct, as we will now show. 

Let us return to the expression (fT5)) for the output as a random variable in the sprin- 
kling process on M 4 . We have already calculated (F) on this basis. There, the three 
random variables in the expression for F are uncorrelated. Now we seek to calculate 
var(F) = (F 2 ) — (F) 2 . The term (F 2 ) can be written 

<F 2 )= g V f fff (a(ds 1 )x(d 4 y 1 )a(ds 2 )x(d 4 y 2 )ax(s 1 ),y 1 )C(x(s 2 ),y 2 )). (32) 

JP JV JP JV (a) 

where si and s 2 are both proper time along P, y\ and y 2 range over the detector D, and 
we have used the symbol (a) as shorthand for the conditions x(si) -< y\, x{s 2 ) -< y 2 . In 
this case, not all of the random variables are uncorrelated. Firstly, x(d 4 y 2 ) is correlated 
to £(x(si), t/i) as z/2 can fall inside the causal interval I(x(s\), y\). In this case if a point is 
sprinkled in d 4 y 2 , <^(x(si), y{) = 0. Therefore there is no contribution to the integral from 
the range in which y 2 € I(x(s\), y\). Using similar reasoning we can put 4 restrictions on 
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the region of integration without affecting its value: 



x{ Sl )<th, (33) 

yiih, (34) 

x{s 2 )£h, (35) 

Viih, (36) 

where we introduce the notation I\ :— I(x(si),yi) and I2 := I{x{s2),y2)- We will use the 
symbol (b) as shorthand for these restrictions. Inside the remaining region of integration, 
all x and o variables are uncorrelated with £ variables. 

Secondly, the x an d o variables can be correlated with each other, as they do not 
all range over disjoint regions. When y\ is in d A y2 there is a correlation. Taking this 
correlation into account is crucial to obtain the correct result (the situation is similar to 
accounting for the self-correlations in a discrete set of variables) . Given a sprinkled point 
at j/i, when 3/2 7^ 2/i ; the probability of finding a sprinkled point in d^y 2 is infinitesimal. 
But when 3/1 — 3/2 the probability is 1, explaining the importance of this correlation. We 
deal with this by splitting the integral into a sum of 4 terms: 

(F 2 ) = Jj + J 2 + J 3 + J 4 . (37) 

These terms take the form 

J t = q 2 K 2 / / / (^(ds 1 )x(d 4 yi)a(ds 2 )x(d 4 y2K(x(si),yi)C(x(s 2 ),y2)) , 

JP JD JP JD (a) (6) (jj) 

(38) 

with the additional (b) bounds from eqns. (|33M36p above, and also some new bounds (ji) 
for each term: 



(ii) 
Ih) 
(h) 
ik) 



s\ = s 2 , yi = 3/2 (39) 

si 7^ s 2 , 3/1 = 2/2 (40) 

si = 52,3/17^2/2 (41) 

si 7^ «2, 3/1 7^ 2/2 • (42) 



In each term, once the obvious identifications have been made there are no remaining 
correlations between x and er variables. Finally there are correlations between the two £ 
variables, which will be dealt with in due course. We will also define 

Ki = J 4 - (F) 2 , (43) 

so that 

var(F) = Ji + J 2 + J 3 + K 4 . (44) 

We now proceed to calculate these four terms. 

Let us start with the J± term. Due to the (ji) bound, which can be used to eliminate 
S2 and 3/2, this term can be written 

Ji=q 2 K 2 [ [ (a(d S ) X (d i y)C(x,y)), (45) 

JP JVx-<y 

= qn(F). (46) 
Consider now the term J2. Calling 3/1 = 3/2 = 2/ we have 

■h = q 2 n 2 f f f (a(ds 1 )a(ds 2 )x(d i yK(x(s 1 ),y)ax(s 2 ),y)) (47) 

JP JP JV (a) (6) 

where the bounds should be interpreted with 3/1=3/2 = y- These bounds are incompatible. 
Bound (a) implies that x(si) ~< y and x(s 2 ) -< y. But all points on P are causally related, 
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including x(s±) and x(si), so x(s 2 ) -< x(si) -< y or x(si) -< x(sa) -< y. This contradicts 
either (f33|) or (f35]l of the (6) bounds. We therefore have 



J 2 = 0. (48) 

The term K4 suffers a similar fate: 
,h = q 2 n 2 f fff (a(ds 1 )x(d i y 1 )a(d S2 ) X (d 4 y 2 )C(x( Sl ),y 1 )<:(x( S 2),y2)) (49) 

JP JV JP JV (a) (b) 
2 .,2 / / / / . ,4 > ,4„ 



= q 2 * 2 / / / ds 1 d 4 y 1 ds2d i y2(ax(s 1 ),y 1 )((x( S 2),y2)) (50) 

JP JV JP JV (a) (b) 

= q 2 K 2 [ fff ds 1 d 4 y 1 d S 2d 4 y2exp{-\I 1 Ul2\), (51) 

JP JV JP JV (a) (b) 

where h is the causal interval between x(si) and yi as before, and |.| indicates the volume 
of a region. The last line comes from seeing that the product ((x(s\), yi)C(x(s 2 ),y 2 ) 
is only 1 when both of these intervals are empty of points, exp(— \I\ U 7 2 |) being the 
probability for this to happen. Subtracting (F) 2 , 

(F) 2 ^q 2 n 2 f f f f ds 1 d i y 1 ds 2 d 4 y 2 exp(-\I 1 \-\I 2 \) (52) 

JP JV JP JV (a) 

gives 

K 4 <q 2 n 2 f fff ds 1 d 4 yids 2 d 4 y 2 \exp{-\hUl 2 \)-exp(-\I 1 \-\I 2 \)\, (53) 

J P JV J P JV(a)(b) 1 J 

where we have added the (6) bound to the (F) 2 integral, which can only increase the 
value of this upper bound on K4. The integrand here is zero unless 1% D I 2 7^ 0- It 
can be seen that there is no region inside the bounds satisfying this condition. As in 
the J2 case, we use the assertion that x(s\) -< x(s 2 ) or x(s2) -< x(s\). Let us assume 
x(s2) -< x(si), without loss of generality, as the conditions are symmetric between Si and 
s 2 . The condition l\ n I2 7^ implies x{s\) -< j/2- Together with x(s 2 ) ~< %(si) this gives 
x(si) £ I 2 , contradicting (f3T)|) in the (a) bounds. Hence, 

K 4 < 0. (54) 

The only remaining term is J3: 

Js = q 2 K 2 f f f dsd^yx d A y 2 exp(-\hUl 2 \), (55) 

J P JV JV (a) (b) 

where now, because of the identification si = s 2 = s, we have I\ = J + (x(s)) D J~(y{) 
and similarly for I 2 . Noting that \I\ U /2I > (|7i| + | J2D/2, we can write 

J 3 <q 2 K 2 f f dsd 4 yi eM~l\Ii\) f A 2 exp(-~|I 2 |). (56) 

JPJV(a) 1 JV(a){b) 1 

Let us define r$ = y| — x 1 (s) and tj — Hi — x°(s), and it.; = fj — fj and Vi — ti + 
(note that in this case we will be holding the a; variable constant while integrating over 
y, conversely to the method by which the mean was calculated). As a step to bounding 
J3, consider 

X(x(s), yi )= f d 4 y 2 exp(-i|/ 2 |). (57) 

This range of integration for 7/2 is much larger than necessary. Let us assume that the 
light-cone of x(s) does not intersect the initial or final spatial boundaries of the detector 
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region (for the points in P at which this does not hold, the region of T> integrated over 
is smaller, so this results in an upper bound on X{x(s),y\) for all points in P). In this 
case, the bounds on the integral are R — x 1 (s) < T2 < R — x 1 (s) + d, < y 2 < \[A and 
< y\ < \f~A. As before we will use T(x(s),y2) 4 ~ (W2U2) 2 . The bound (a) gives U2 > 
in this case, and so we have 



r 2(R-x 1 +d)+u 2 
>2(R-x 1 )+U2 

As before we have used r(s, y2) i ~ {U2V2) 2 ■ We have 



A POO r2(R-X 1 +d)+U 2 

X(x{s), yi )<- du 2 dv2exp(- — (u2V2) 2 ). (58) 

1 JO J2(R-x 1 )+U2 48 



X(x(s), yi )<Adj du 2 exp(-^(u2{2(R-x 1 ) + U2)) 2 ) (59) 

< VSAd]-, (60) 
R 

since d and x 1 are negligible compared to R and as in the calculation of (|!2T|) we have 
neglected insignificant terms in the exponential. Inserting this back into (|56[) we have 

J 3 <V3q 2 K 2 Ad±- [ [ dsd 4 y 1 exp(-i|/ 1 |). (61) 

K JPJV(a) 1 



This integral is similar to the one we have already calculated to find (F) in eqns. 
We find 

J 3 <VeqKAd-{F). (62) 
R 

Finally, summing up the remaining terms in eqns. (|46|) and (|62[) to find var(F) from 
eqn. (|33)) : 

var{F) < q(F)(l + V^nAd-). (63) 

R 

Now, var(F) represents the variance of the signal from only our oscillating positive 
charge. We are not actually interested in var(F) but rather in the variance of the total 
detector response to the oscillating positive charge andjhe static negative charge together, 
far(Ftotai) = var(F + +F-), where we have renamed F above as F + as before. This total 
variance is 

var(F tota x) = var(F+) + var(F_) + 2 Cov(F + ,F_), (64) 
where Cov is the covariance function, 

Cov{F +1 F_) = ((F + - (F+))(F- - (F-))). (65) 

It is a standard corollary of the Cauchy-Schwarz inequality that 

Cov(F+,F^) 2 < var(F+)var(F-), (66) 

and so we have 

var(F tota i) < var(F+) + var(F-) + 2\J var(F+)var(F- ) (67) 

< ivar(F + ) 7 (68) 

where the result follows from (|63|) and (|29|) : var(F-) is the same function as var(F+) 
with the source velocity v set to zero so var(F-) ~ var(F+) for low velocities. 
The standard deviation of the total signal as a proportion of the output is 

\var(F totai ) \Uvar(F + ) 

-i — = 1—^= . (69) 

(■Ptotal) (-Ptotal) 
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Using our result ([63]) for var(F + ) (and the value set earlier for k) gives 



var(F to t a i) 

(^total) 2 



< 16 



v 



■ 2 TAd 



ttR 



+ 



2v 2 T 



1 



(70) 



We can now input some values for the parameters to estimate the fluctuations for a 
model source and a model detector. We will obtain an upper bound on the fluctuations by 
making unrealistically restrictive estimates in some cases. For a description of a relevant 
source, we look to discussion of coherence of electromagnetic radiation from extragalactic 
sources in the literature. One coherent source that has been cited [10] is the Active 
Galactic Nucleus PKS 1413+135 [23]. This source is at a distance of order one Gigaparsec, 
and coherence was detected in radiation of wavelength 1.6 /im. Therefore R is set to be 
of order 10 60 Planck lengths, and in order that the time resolution of detection be much 
better than the period of the radiation (one of the unrealistically restrictive estimates) , we 
set T to be 10 27 Planck units. The radiation from AGNs is thought to be synchrotronic, 
meaning that the source electrons have velocities close to 1. This doesn't fit with our low 
velocity assumption but we choose ^0.1 as a compromise value for which our expansions 
are still approximately valid d . As for the spatial size of the detector, y/A and d can 
safely be set to be of order 1 nanometer, or 10 26 Planck units. With these values of the 
parameters the first term inside the brackets in equation (|70|) is of order 10 -43 and the 
second is of order 10~ 25 . With these values the standard deviation of the total signal is 
about one part in 10 12 . This order of fluctuation is well below that which is possible to 
detect even in more realistic set-ups than this nano-size and femto-duration detection. 

This shows that, with these values for the parameters of our model, the fluctuations 
are a small part of the total signal. Thus for a typical sprinkling, it is justified to claim 
that the results of the discrete model match the continuum model sufficiently well for the 
difference to be undetectable. The conclusion is that the spatial and temporal coherence 
of waves from distant sources is preserved in this model. 

It is interesting to ask if models of this type can give good results for all types of 
source. In this connection, it should be noted that certain approximations and assump- 
tions, mostly unfavourable to the eventual outcome, have been made to achieve the result 
given above for the distant AGN. The bound on fluctuations coming from equation ()70D 
is well above that which could reasonably be achieved in such models, for three main 
reasons. Firstly, the values for the quantities involved are all set to give a higher result 
for the fluctuations than is realistic, to give an upper bound; the only one that might 
conceivably be less favourable for other sources is v. Secondly the bound on var(F + ) 
given in (l63|) is fairly crude: not only could a smaller upper bound be found for J3, 
but the K4 term could be less than zero. Thirdly, making the model more realistic, by 
considering detection lasting for multiple oscillation periods, and/or spatially extended 
source^, would be expected to greatly reduce the relative size of the fluctuations. Thus it 

6 To fully satisfy scepticism on this point, we would have to repeat our calculations without this 
approximation. This would necessitate dropping the further assumption that the duration of the detection 
is much shorter than the period of oscillation of the source: when the maximum source velocity is close 
to the speed of light, the field from the postive charge would vary from much greater than its value for 
v = to much smaller, and so it is possible that the fluctuations as a proportion of the dectector output 
would vary greatly over one oscillation. Such a large value for the relative fluctuations at only one time 
in the oscillation would not in practice be measureable, however; it is more realistic to consider a longer 
detection for the cited source in any case. Considering this, there is strong reason to think that this 
more complicated treatment would not lead to a different conclusion. When v is positive, the ratio of 
fluctuations to detector output become even smaller as v increases towards the relativistic regime, and 
when vel w — 1 the contribution to the total detector output (and fluctuations) will be small. 

7 An extended source could be built by adding many positive and negative particles, in which case a 
bound on the total variance similar to l|67jl would still hold. The resulting bound would be proportional 
the square of the number of particles, and -Ftotal would of course be proportional to the number of 
particles, leaving H70H unchanged. However, in this case the result analogous to would be a very 
much looser bound, as signals generated from significantly separated source particles would not be highly 
correlated. The relative fluctuations in such a model would arguably be greatly reduced as a result. 
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is reasonable to conjecture that models of this type could be constructed to model many 
realistic electromagnetic sources^]. In any case, the main aim here was to show that the 
model is consistent with the coherence of waves from distant sources. 

4 Conclusion 

We have reviewed a simple model of the propagation of a scalar field from source to de- 
tector in a Minkowski spacetime, and demonstrated a model of the same physical system 
in which the continuous spacetime is replaced with a causal set. The main result is that 
the output of the detector is not significantly affected by the introduction of discreteness. 
This being so, the phase of the source as it oscillates is faithfully recorded at all astro- 
physical and cosmological distances: there is no loss of (classical) coherence. Nor is there 
any change in the dispersion relations. The discreteness will eventually show up when the 
distance, i?, between source and detector becomes large enough that the mean number of 
links from source to detector region is of order one but the signature of the discreteness 
will be a "cutting out" of the signal altogether rather than any loss of coherence. Un- 
fortunately, for this to happen for any reasonably sized detector, R would have to be a 
super-horizon distance and in any case the signal would by then be undetectably tiny. 

The relevance for quantum gravity phenomenology of a model premised on a fixed 
background structure may be questioned. After all, the continuum seems likely to re- 
sult from a "quantum superposition" of causal sets, whatever that turns out to mean. 
However, the purpose of this article was to show that discreteness per se does not pose 
any problem for the coherent propagation of waves over long distances. Claims about 
"generic" phenomenological predictions from quantum gravity will need to take account 
of this. For example, it cannot now be claimed that a fuzziness in spacetime metric 
relations necessarily disrupts Huygen's principle nor that modified dispersion relations 
always result from the existence of a minimum length scale. We emphasise, again, that 
the maintenance of coherence and the usual dispersion relations hold for any fixed causal 
set that arises with relatively high probability from the sprinkling process and is not the 
result of any averaging over causets. 

The crucial aspect of causal set discreteness that distinguishes it from other approaches 
is that it respects Lorentz invariance. This is the key to preserving the usual dispersion 
relation in the above model: even without a detailed model, we would have claimed that 
variance in the speed of light as in 10 could not be present in a causal set model unless 
it were introduced "by hand" somehow. The concrete model is particularly illuminating 
because it shows explicitly that the speed of photons from distant sources is not affected 
by causet discreteness due to the existence in the causet of structure - the links - which 
traces lightcones extremely accurately. It is hoped that the study of Gamma ray pulses 
from Blazars and Gamma ray bursts [HE] will either detect, or place stringent limits on, 
variations in c and it remains to be seen which discrete models of propagation will be 
consistent with these observations. 

We have thus seen that, due to the tininess of the Planck scale and the Lorentz in- 
variance of the approach, causal set models can give results that are experimentally indis- 
tinguishable from the continuum. This is simultaneously encouraging and disappointing. 
Encouraging because every such result provides more evidence for the "Hauptvermutung" 
(Central Conjecture) of causal set theory, that causal sets can indeed be the deep struc- 
ture of Lorentzian spacetimes. The disappointment arises because providing experimental 

8 Indeed, the situation might reasonably be expected to be better for a similar discretisation of the 
electromagnetic Green's functions, if the calculations were otherwise similar to those above. In the 
continuum, for the scalar field, the DC and AC parts of the signal have the same order of magnitude 
when v ~ 1, whereas in the electromagnetic case the the AC part will be larger by a factor of R. This 
would suppress the total variation of the signal but not the total mean signal. This could be analysed by 
investigating the causet version of the continuum model suggested at the end of section 2 in which the 
field 4> is treated as an electostatic potential. 
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evidence for causal set theory will require predictions that are at variance with those made 
in the continuum. Though the simple model presented here seems not to be able to pro- 
duce such phenomenology, other models of matter propagating on causal sets do hold 
out the hope of producing observable signatures of discreteness. One promising class of 
models is based on the idea that if spacetime is a causal set then this should cause small, 
random, Lorentz invariant fluctuations in the motion of particles through spacetime. For 
massive particles this results in a diffusion in momentum [5j [25] and for massless parti- 
cles, both a diffusion and drift in frequency [25]. These effects are potentially observable 
and astrophysical and cosmological data has been used to bound the parameters of the 
models [5J US] [371 Beyond these models, the discovery of a discrete, Lorentz invari- 
ant, D'Alembertian operator p~3j [14] means that a discrete scalar wave equation can be 
solved and the behaviour of waves and wave packets investigated, providing more realistic 
models of matter and further opportunities for predicting phenomena that could reveal a 
fundamental spacetime discreteness. 
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